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ABSTRACT 



AN INVESTIGATION OF THE VIRTUAL MASS OF A CYLINDER 
VIBRATING IN WATER 

bj 

DAVID A. ROGERS MCLEAN C. 3HAKSE0BER 

LIEUTENANT, U.S. NAVY aoA LIEUTENANT, U.S. NAVY 

Submitted to the Department of Naval Architecture and Marine Engineering on 
May 25» 1953 in partial fulfillment of the requirements for the degree of 
Naval Engineer. 



The object of this thesis is to experimentally investigate the virtual 
maos of a hollow cylinder vibrating in water. 

A lucite cylinder was magnetically vibrated in air and water at various 
length to diameter ratios and the frequency of vibration for as many modes 
as possible, up to five, recorded. No attempt was made to measure ampli- 
tude. The ratio of added water mass to displaced water mass was computed 
from the frequencies and compared with analytical results. 

The investigation shows that end effects have a very great influence 
on virtual mass. A9 the length to diameter ratio is decreased, the added 
virtual water mass is decreased. There is also a decrease in virtual water 
mass as frequency is increased at constant length to diameter ratios. 

A ratio of the measured virtual water mass to analytical, called K, was 
computed and found to be a function of L/D and mode number. 

It Is recommended that further investigations using bodies of revolution 
whose ends have zero area such as ellipsoids be made. It would be desirable 
to use equipment to permit measuring amplitude as well as frequency. 

Where data is available, it is recommended that an attempt to calculate 
the frequencies of an actual hull such as a submarine be made, correcting the 
analytical virtual water mass by the applicable K values. 



Thesis Supervisors Professor F. M. Lewis 
Title* Professor of Marine Engineering 
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it introduction 

l 

When laser sed in a dense fluid, such as vater, a body vibrates as 

though it had undergone an increase in mass. This increase in mass is due 

to the flow of fluid about the body as it moves* 

Several analytical methods of computing the virtual mass of a body 

(1 2 ) 

vibrating in a dense fluid have been presented*' * ' These methods are 

based on the assumption of potential flow about the body* Unless the body 
is of uniform shape with pointed ends, i.e. an ellipsoid, it is not presently 
possible to compute analytically the virtual mass because of the flow about 
the ends being highly rotational. 

A submerged submarine must vibrate as a free -free body at a frequency 
determined by its virtual mass as described above. When vibrating horizon- 
tally, its virtual mass will be somewhat lower than that computed on the basis 
of no end losses because of flow about the ends. 

Professor Frftnk M, Lewis has presented a method for determining the vir- 
tual mass based on an ellipsoid which may be corrected for other shapes. ^ 

Ur. H. M. Schauer of the Underwater Explosion Research Division, Norfolk 
Naval Shipyard, has done likewise for a oylinder with no end flow. Mr. E. B. 

Koullin and Kr. A. D, Browne, in a paper presented before the Cambridge 

(3) 

Philosophical Society in 1928 ' gave the results of their investigation of 

the periods of a free-free bar of rectangular cross section vibrating in water. 
In their experiments they used long bars which had length to depth ratios of 
from 26 to 39. Using such long bars they found that flow about the ends did 
not have any appreciable effect on the virtual mass as analytically computed. 
However, they did not investigate lower length to depth ratios. They found 
that the virtual mass is not affected by depth when below about two diameters. 
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The following Is a report on the experimental determination of the 
virtual mass of a circular cylinder for several length to diameter ratios 
while vibrating in water* 
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II. PROCEDURE 

In this section the steps followed to obtain the desired results are 
described. The procedure consists of two parts, experimental and analytical* 
Details of the procedure arc presented in the Apoendix. 

Experimental Procedure 

Proa "The Theory of Sound" by Rayleigh^ the appropriate equations 
were used to obtain the nodal and anti-nodal points of a free-free bar. 

Using the frequency equation for a free-free bar 




the first five modes were computed to give an approximation of the natural 
frequencies. 

To vibrate the cylinder mechanically would require a motor with a speed 
range of 1800 to 60,000 revolutions per minute. For this reason, magnetic 
vibration of the bar was by far the preferred method. Schematics of the 
apparatus are shown in Figure I. 

A Lucite plastic tube 52.7 inches long with an outer diameter of 2 inches 
and an inner diameter of 1.75 inches was used for ths first test. On the 
end of the cylinder was wrapped same small diameter soft iron wire to permit 
magnetic excitation. The amount of wire was not great enough to affect the 
frequency or mass of the bar. By experimenting with various types of pickups, 
it was found that a seismic crystal gave the best results* The pielcup was 
very light in weight and very sensitive to vibration. This particular pickup 
was a Brush seismic crystal used on the sounding board of an electric guitar. 
The pickup was mounted on the inside of the cylinder at the opposite end of 
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the exciting wire. A small hols was drilled at a nodal point through which 
the vires freen the pickup were run* 

An audio osoillator with a frequency range of 20 to 20,000 cycles per 
second vac used to drive an electro-magnet* The output of the audio oscill- 
ator was amplified in a power amplifier and this was used to drive the elec- 
tro-magnet* The output of the crystal pickup was put into a cathode ray 
oscilloscope where the signal was peaked for resonance* Unfortunately, due 
to radiation, frequencies above 1600 cycles per second could not be detected* 
Two different types of frequency measurement were employed to compute 
the frequency of vibration, both of which gave very accurate and like results* 
The first method was to take the output of the pickup and put it on the verti- 
cal plates of the CRO* When a resonant signal was obtained, a picture of the 
frequency was taken with no horizontal sweep* The camera used was a very 
high speed model with no shutter* The camera had a built-in timing light 
of 1000 cycles per second vhioh shoved on the film. The frequency was then 
computed from the developed film by counting the cycles* 

The other method of determining the frequency was to take the amplified 
output of the crystal pickup and put it Into an electronic decade counter* 

To compute the air frequencies, the cylinder was suspended by strings 
located at the nodal points* The electro-aagnet was placed as close as poss- 
ible to the soft iron windings. To prevent banging of the cylinder, a rubber 
band was used as a standoff* The frequencies were recorded as described 
above. 

For the water teats, the cylinder was immersed seven diameters in the 
towing tank in the M» I* T. Hydrodynamics Laboratory. This depth ensured that 
no surface effects would be present. The cylinder was anchored by two strings 



at the nodes 
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The sane procedure was followed U8lng cylinders of 40, 33, 30, 23, and 



22 inches, and for a 42-inch bar with six-lnoh conical ends. 



The ratio of the added water aass to displaced water mss, M yA*- Q i was 
ooRputed directly fro® the observed frequencies as explained in Details of 
Procedure. 

A correction factor K wa» then calculated, where K is the ratio of 

(o) 

neasured M to that ccouuted by r r. H. H* Sehauer. ' 1 
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III. RESULTS 



The results are plotted as shown in Figs* II through V* 
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IV. DISCnSSIOK OF RESULTS 

The results shown plotted In figures II through V indicate that the 
ends have a largo influence on the value of the ratio of added water mass 
to the displaced water mass, i,e. M As the length to diameter ratio is 
decreased, the value of K, (the ratio of measured M to the M w /^ c cal- 
culated by Dr. Schauer) , decreases rapidly for a given value of f|* However, 
K Increases for increasing f^&t constant L/D» 

K is based on an analytical curve which assumes a potential flow, and 
therefore irrotational, with no flow about the ends. Thus K may be con- 
sidered a measure of the end flow and of rotationality. For a cylinder with 
flat ends a large portion of K can probably be attributed to flow of the 
fluid about the ends. 

Just how much of the difference from analytical values can be attrib- 
uted to end flow is difficult to say from these experiments. An inspection 
of the test results of the cylinder with conical ends indicates that end 
flow is not the only contributing factor since the area at the ends of the 
cones is substantially zero. However, due to the abrupt change in shape of 
the body where the cones are Joined to the cylinder, a certain amount of 
spilling of fluid toward the ends will occur. 

Due to the large effects of L/D on the value of M it appears that 
a direct application of Dr, Schauer' s equation is not feasible unless the 
body is very long relative to ita breadth. 

Professor F. M» Lewis ^ has shown that the added water mass per 



foot of length for an ellipsoid is 
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m • cjtr» 2 jr 

w * w 



where C * Section Inertia coefficient 

B * Half bean or radius of ellipsoid at section 

« Specific weight of fluid 

- m Actual K. E» surrounding fluid _ 

K. E. of fluid if flow Is two dimensional 

Now for a cylinder C » 1 and we can write 





H * JtTE y L 

V M V 


(3) 


and 


M JirR 2 * L 

Jf * J-iL 

© If 


(4) 

' i 


so that 


M 

J *> ^ 

c 


(5) 



It was found that if Prof, Lewis* J factor is multiplied by our K, a 
fairly decent approximation to the actual frequency of a submerged body can 
be computed, A computation was made of the submerged frequency of a 42-inch 
body having 6-ineh cones at each end of a cylindrical middle body 30 inches 
long, 

Using Prof, Lewis’ J factor multiplied by K for an L/D ratio of 21 at 
= ,15* the computed two noded frequency was 68,9 c.p»s» and the experiment- 
ally measured value was 70 e,p.e» For the three noded frequency the computed 
value was 187 c,p.s. and the measured value was 188 c,p.e» 

Although this method gives very good results in the case tested, further 
investigation is required to determine its limits of application. Similar 
results would have been found using Dr* Schauer ’ s in a like manner* 
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V. KEC0?fl<EM)ATI0B3 

In view of the large discrepancy between values of theoretical mass 
found by analytical means and those measured. It is recommended that an 
investigation similar to this be made of families of bodies of revolutions 
having ends whose section decreases to zero. 

The measurement of the virtual mass of bodies which do not have complete 
radial symmetry is much more difficult since pains must be taken to insure 
that the vibrations are limited to the plane desired. However, it is 
recommended that where data exists on the vibration of submerged bodies such 
as submarines, an attempt be made to apply the correction, K, to compute 
their frequencies. 
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Selection of Tubing 

Sine® the object of this investigation is to determine the virtual 
added mass of water, it is desirable that the cylinder be made of a material 
having a low specific gravity. Thus the added water maos will be & large 
fraction of the total virtual Bass of the body. In order that the frequen- 
cies of the cylinder will not be excessively high, the material should have 
a relatively low modulus of elasticity* 

As a result of these considerations, it was decided that the cylinder 

should be constructed of Incite. The specific gravity of Luc its is 1.18 

5 

and has a modulus of elasticity in the order of 5 x 10 pounds per square 
inch. 

Effect of Added Hasses 

The weight of the soft iron wire was .373 ounces and covered three- 
fourths of an inch of the cylinder. The weight of the crystal pick-up and 
clamp used to hold It securely In place wae .303 ounces. The weight of the 
pl&ntic cylinder was .497 ounces per inch* The effect of these added masses 
on the frequencies was assumed negligible. 

The high speed movie camera method of frequency measurement proved very 
satisfactory, but also required a great deal of time. The camera had a neon 
bulb timing light built in which was energised by a 1000 cycle audio oscilla- 
tor. The output of the oscillator was amplified through a CR0 which provided 
a means of adjusting the light to the d*3ired brightness. The only errors 
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that could be made in this method wore (a) error# in counting the film and 
(b) frequency drift in the oscillator* Two separate runs wore made on each 
length tested. 

In order to speed up the experiment and reduce the amount of labor 
involved, it was decided to try frequency measurement with an electronic 
decade counter* If there ia any frequency drift of the audio oscillator, 
the results of the decade counter will be more accurate than the movie cam- 
era since the decade counter records every second noting any change which 
may occur. 

Boundary Eff ects in Water 

To check for wall effects, the cylinder was submerged in the stability 
tank and the test rerun. The results obtained war® the same ae those 
observed in the towing tank teat* A further check vae made by varying the 
distance of the cylinder from the walla in the stability tank. Again no 
difference was noted. Ste attempt was made to vibrate the cylinder within 
four diameters of the wall* 

Comparison with Different Cylinder 

The necessity of re-checking the 30-inch cylinder arose after it had 
been cut down to 22 inches* A new 30-inch plastic tube was re-wound with 
approximately th# same amount of wire as before* Both the air and water 
tests agreed identically with the previous 30-inch test. 

Damping Jffects of Bn bber Band Stand-off 

To ascertain whether the rubber bard had any appreciable damping effect, 
the force necessary to stretch the rubber band one inch and the force required 
to deflect the tube one inch were calculated* To stretch the rubber band 
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on# inch required n force of *61 pounds, and to defleot the cylinder one 
Inch, a farce of 454 pound a was n«o« 9 nery. The damping effect of the rubber 
bend arr be neglected* 

feaUMff 2ln&m:£*c 

The equipment used wee not capable of sceaeuring amplitudes® Thie ie 
partly because the natural reaoajmt period ef the pick-up, 980 cycle® per 
second, was within the range of the frequencies measured® The primary 
difficulty was that it was not possible to maintain a constant oxciting 
force on the cylinder* There also was no roan* of accurately measuring the 
exciting force* A constant exciting force was not possible because a higher 
force was required to excite the higher modes $ but this same force would 
cause the cylinder to he drawn hard against the faces of the electro-magnet 
pedes at tho lower nodes* It was necessary therefore to start with relatively 
low driving forces at low nodes cad increase the force to excite tho higher 
frequencies* 



The virtual mss was computed assualng negligible damping so that 
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where 

• Virtual Ms.sc 
Hj, *• V&nn of Tube 

? 1 * *» *ir 

fg * Freqtiency in Water 
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M q ■ Displaced Water Ma 3 s 
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Thus the ratio ( ) may be computed from the measured values of the 

% 

frequencies* 



( 7 ) 

( 8 ) 

( 9 ) 



Computation of K 

( 2 ) 

Dr. H. M. Schauer has derived an analytical expression for the mass 
M 

ratio, , as follows* 
c 




where H Q and are Hankel functions and other symbols as previously defined* 
The variation of the measured mass ratio from the above will be a func- 
tion of 1^ and the ratio* This ratio can be expressed as a ratio, K 

5 „ measured mas 3 ratio 
analytical mass ratio 
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B » SPtMABI OF DATA AND CALCOTATIOfO 

1. Definition of Symbols. 

f^ * Observed air frequency, cycles per second 
» Observed water frequency, cycles per second 
■ Added mass of water 
* Mass of displaced water 
L » Length of cylinder, inches 
D » Diameter of cylinder, inches 

>1 

m » Mode number 
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2. Calculated 'rallies of ^ 

c 
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f l 


f 2 


M 

c 
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1 


61.7 
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2 


170.5 


86.2 


0.82 


0.179 


3 


332.0 


169.1 


0.795 


0.238 


4 


539.5 


279.0 


0,770 


0.298 


5 


787.5 


419.0 


0,709 


0.358 



TABLE I 

Teat Results 52.7-inch Cylinder L/b * 26.35 
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.682 


.471 



TABLE II 

Test Results of 40-inch cylinder L/b = 20 
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ItCDE 


f l 


f 2. 


M . 


1 


1 


114 


59 


0,773 


0.166 


2 


313 


163 


0.761 


0.249 


3 
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319 


0,727 


0.332 


4 


943 


509 


0.665 


0.415 


5 


1351 


739.5 


0,665 


0.497 



Test Results 3&-inch Cylinder L/D * 19 
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0.523 
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MODE 


f i 


*2 


M o 


1 


1 


205 


112 


0*674 


0.224 


2 


553 


305 


0*657 


0.336 


3 


1036 


551.5 


0.622 


0.448 
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Test Raeults 23-ineh Cylinder L/D * 14 
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Teat Results 22-lnch Cylinder I./b « 11 
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*120 


1.431 
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,275 
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1.0330 


,968 
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Analytical Calculation of Virtual Hass 



By Dr* H» M, Schauer ^ 




H ■ Added Water Haas 
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1 






« (1^) 


tyij) 


iM_(ii) 

1+ 1 


M 

K 

0 


,120 


1.431 


5.20 


,275 


,0330 


1.0330 


,968 


131 


1.376 


4.75 


,290 


.0380 


1,0380 


.963 


175 


1.193 


3.50 
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.0599 


1.0599 


.943 
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Analytical Calculation of Virtual Mass 
By Dr» H, M. Schauer^ 




■ Added Water Maas 

■ Mass of Displaced Water 
fij, » Macs of Cylinder 
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a * Mode Number 
a » Radius 



L * Length 
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3. Calculated values of K 
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<Wh 


(h A>a 


K 


■ 2 


.811 


.933 


.870 


.25 


.792 


.910 


,870 
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.772 


.886 


.870 


35 
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.860 


.870 
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.833 


.876 






TABLE IX (a) 






Values 


of K L/b « 26.35 
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Values of K L/D * 19 
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1 
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L/D - 15 
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TAPLE IX (e) 
Values of K L/D * 11 
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Values of K L/fa * 21 
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C. SAMPLE CALCULATIONS 



1. Theoretical frequency of free-free bar vibrating in air. 



f 





1/2 



where 



* * 4<*V30 for first mode 

21 

E • Modulus of elasticity 

• 580,000 psi for Lucite (appro*.) 
Y * Specific gravity * 1.18 
g «* 386 in/see 2 

k * Radius of gyration of section 

k 2 » .441 



:QQ81 ^.QSQJS-tAfl 
6.28 11.05 x io-5 



1/2 



* 62*0 cycles per second 
for 2 noded frequency of 52. 7-inch cylinder. 



2. Calculation of 
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■ Added water mass 

* Displaced water mass 
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For second node of 52.7-inch bar 

Wy * Wt/in. x L + Wt. clamp + Wt, Pickup ♦ Wt» Exciting wire 

u , 

16 x 32,2 

* 0.0521 

M 

C 4g 

4 S 32,2 x 1728 

* 0,186 
f x « 170,5 
. f 2 * 86.2 
M 

r* - , 280 ( 2 , 92 ) 

c 

- 0,820 

3. Calculation of K, 

By definition 

M v/ * c measured 

^ ** M^/M~ analytically computed 

For 52.7" cylinder 
when n 3 .2 

(^)H - * Bli 

0 

<^>A * - 933 

Heno. K » - 0,870 
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The center of gravity of cones is approximately 3" fro® the base. Let 
ns assume this body to be equivalent to a right circular cylinder 36 inches 
in length* 

This assumption may be checked by calculating the air frequency of a 
36»inch cylinder. 



where 



r A ( kx ) 1/2 

21T ' y * 



mr2_ 1/2 

( “jr® ) * 7,350 



2ir ' * 

*a * f or 2 noded frequency 

f - (Hjf 2 ) 2 7350 

f * 127 ops 

Observed frequency was 131 ops 9 therefore equivalent length is 



35»5 inches. 
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The w ter frequency now may be computed using Prof. Lewis* J, corrected 
by K for a 42-inch right circular cylinder 



f . 



f ( 7350 ) ( 



,/z 

V H v 



where K^, *= ,036 

M « JKM 

V c 

« (,947) (,839) (,120) 
\ * *0952 

h w - *1312 



f * 68,9 eps. 

The experimental results gave a water frequency of 70 cpa. 



( 12 ) 



For the second mode 



J * • 923 

K « ,849 

K y « ,0940 

M y ♦ Jtp * .1300 

f = 187 cpa 



Experimental results gave a water frequency of 188 cpa 
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D. ORIGINAL DAT/ 



1* The following data was taken using the decade counter. 



MODE FREQUENCIES AVERAGE 



1 


114 


114 


114 


114 


2 


313 


313 


313 


313 


3 


604 


601 


602 


602 


4 


944 


942 


942 


943 


5 


1352 


1350 


1351 


1351 



Air Frequencies 38-inch Cylinder 



MODE 




FREQUENCIES 




AVERAGE 


1 


59 


59 


59 


59 


2 


163 


163 


163 


163 


3 


320 


319 


319 


319 


4 


510 


507 


509 


509 


5 


740 


739 


739o5 


739.5 






TABLE XI 







Water Frequencies 38-inch Cylinder 
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MODE FREQUENCIES AVERAGE 



1 


182 


183 


183 


183 


183 


2 


496.5 


495.5 


496.5 


496.5 


496.5 


3 


933 


936 


935*5 


936.5 


936.0 


4 


1454*5 


1453*5 


1453 


1454*5 


1454.0 



TABLE HI 

Air Frequencies 30-inch Cylinder 



MODE 




FREQUENCIES 




AVERAGE 


1 


98 


97.5 


97.5 


97.5 


97.5 


2 


268,5 


269 


267 


268,5 


268.5 


3 


510 


513 


513.5 


511 


512.0 


4 


810 


809 


808 


807.5 


808 



TABLE HI! 

Water Frequencies 30-inch Cylinder 
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MODE 



FREQUENCIES 



AVERAGE 



1 


205 


205 


205 


2 


553 


552.5 


553 


3 


1036 


1035.5 


1036 



TABLE XIV 

Air Frequences 28-inch Cylinder 



MODE FREQUENCIES AVERAGE 



1 


112 


112 


112 


2 


305 


305 


305 


3 


581.5 


581.5 


581.5 


4 


916 


917 


916.5 



TABLE XV 

Water Frequencies 28-inch Cylinder 
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MODE FREQUENCIES AVERAGE 



1 


298 


299,5 


298,5 


299 


2 


800 


801 


800 


800 


3 


14 57 


1453 


1457 


1457 
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Air Frequencies 22-inch Cylinder 



MODE FREQUENCIES AVERAGE 



1 


171 


171 


172 


171 


2 


465 


466 


465 


465 


3 


854 


355 


855 


855 



ms&jsm. 

Water Frequencies 22-inch Cylinder 
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f'-ODE FREQDENGIES AVERAGE 



1 


104 


105 


105 


104 


104.5 


2 


291 


290 


291 


290 


290.5 


3 


561.5 


562 


562 


562 


562 


4 


891 


894 


893 


893 


893 


5 


1278 


1276 


1276 


1277 


1277 



TABLOmi 

Air Frequencies 40-inch Cylinder 



MODE FREQUENCIES AVERAGE 



1 


54.0 


54*0 


53.5 


54*0 


54.0 


2 


150.0 


151.0 


151 


151 


151 


3 


296 


294 


295 


295 


295 


4 


477 


476,5 


476 


476 


476 


5 


693 


691 


692 


691.5 


692 



TABLE XIX 

Water Frequencies 40-inch Cylinder 
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HODS 




FREQUENCIES 




AVERAGE 


1 


131 


131 


131 


131 


2 


353 


351 


352 


352 


3 


659 


658 


658 


658 


4 


1016 


1020 


1014 


1017 


5 


1417 


1419 


1415 


1417 













Air Freqr*encid3 30-inch Cylinder 
with 6-inch Conical End a 



MODE FREQUENCIES AVERAGE 



1 


70 


70 


70 


70 


2 


188 


188 


188 


188 


3 


358 


355 


352 


355 


4 


555 


555 


555 


555 


5 


796 


790 


790 


792 



Water Frequencies 30-inch Cylinder 
with 6-Inch Conical Eads 
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2* The following data was taken using the high speed movie camera* 



MODE FREQUENCIES AVERAGE 



1 


61.7 


61.7 


61.7 


2 


170.5 


170.5 


170.5 


3 


332.0 


332.0 


332.0 


4 


539.5 


539.5 


539.5 


5 


787.5 


787.5 


787.5 



TABLE mi 

Air Frequencies 52.7-inch Cylinder 



MODE FREQUENCIES AVERAGE 



1 



2 


86,2 


86*2 


86.2 


3 


169.1 


169*1 


169.1 


4 


279.0 


279,0 


279.0 


5 


419*0 


419.0 


419.0 



TABLE mil 

Wqter Frequencies 52* 7-inch Cylinder 
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TWO NODEP WATER FREQUENCY 




TWO MOPED AIR FREQUENCY ' 
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FIGHRE VII 

PHOTOGRAPHS OF TYPICAL SIX NCDED FREQUENCIES 




SIX NOPED AIR FREQUENCY 




SIX NODED WATER FREQUENCY 
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